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<^ ' 1. Introduction 

Cp ■ The discretization of integrable systems is an important issue in mathematical physics. 

The most common situation is that in which some or all of the independent variables are 
discretized. A discretization process in which the dependent variables are also discretized 
in addition to the independent variables is known as "ultra-discretization". One of the 
most important ultra-discrete soliton systems is the so-called "soliton cellular automaton" 
(SCA) [12, 16, 17]. A general method to obtain the SCA from discrete soliton equations was 
K^ ■ proposed in Refs. [6, 18] and involves using an appropriate limiting procedure. Another 

issue which has received renewed interest in recent years is the phenomenon of soliton 

''^ • resonance, which was first discovered for the Kadomtsev-Petviashvili (KP) equation [8] (see 

C^ \ also Refs. [7, 11]). More general resonant solutions possessing a web-like structure have 

recently been observed in a coupled KP (cKP) system [4, 5] and for the KP equation itself [1]. 
In particular, the Wronskian formalism was used in Ref. [1] to classify a class of resonant 
solutions of KP which also satisfy the Toda lattice hierarchy. It was also conjectured in Ref. [1] 
that resonance and web structure are not limited to KP and cKP, but rather they are a generic 
feature of integrable systems whose solutions can be expressed in terms of Wronskians. 

The aim of this paper is to study soliton resonance and web structure in discrete soliton 
systems. In particular, by studying a class of soliton solutions of the two-dimensional Toda 
lattice (2DTL) equation, of its fully discrete version, and of their ultra-discrete analogue which 
was recently introduced by Nagai et al [9, 10], we show that an analogue to the class of 

f E-mail: maruno@math.kyushu-u.ac.jp. 
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solutions studied in Ref. [1] can be defined for all three of these systems, and that a similar 
type of resonant solutions with web-like structure is produced as a result in all three of these 
systems. To our knowledge, this is the first time that resonant behavior and web structure 
are observed in discrete soliton systems. These results also confirm that soliton resonance and 
web-like structure are general features of two-dimensional integrable systems whose solutions 
can be expressed via the determinant formalism. 

2. The two-dimensional Toda lattice equation 



We start by considering the two-dimensional Toda lattice (2DTL) equation. 



dxdt 



K+i - 2Vn + K-1 , 



with Qn{x, t) = log[l + Vn{x, t)]. Equation (2.1) can be written in bilinear form 



through the transformation 

Vn{x,t)-- 



dxdt 



logTn{x,t). 



(2.1) 



(2.2) 



(2.3) 



It is well-known that some solutions of the 2DTL equation can be written via the Casorati 
determinant form Tn — Tn [3], with 



r(N) 



f, 
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(2.4) 



In+N-l 

where {/„ '{x,t), . . . , fn {x, t)} is a set of N linearly independent solutions of the Unear 
equations 



a/, 



(i) 



= ,/, 



df, 



(0 



n+l : 



= -fn- 



(2.5) 



dx ■"'^^' dt 

for 1 < z < iV. [Note that the superscript "(^)" does not denote differentiation here.] For 
example, a two-soliton solution of the 2DTL is obtained by the set {/'^^\ f'-^^ }, with 



/(^)^c^^'-%e«^\ 



i = 1,2. 



(2.6) 



J = l,...,4, 



(2.7) 



where the phases 0^^^ are given by linear functions of {n, x, t): 

el^'^ {x, t)^n \ogpj + pjx — 1 + e[^^ , 

Pj 

with pi < p2 < Ps < Pi- Equation (2.6) can be extended to the A^-soliton solution by 
considering {/^^', • • • , f^'^'}, with each f^^ defined according to Eq. (2.6). 

On the other hand, solutions of the 2DTL equation can also be obtained by the set of 
' '^^ — 1, . . . , M} with the choice of /-functions. 



r functions {r„ ' '' | N 



/, 



ii) 



fn+i-l 



l<i< N < M , 



with 



M 



(2.8) 
(2.9) 
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and with the phases dn , 1 < i < M still given by Eq. (2.7). (Note that the meaning of N 
and M is the opposite of Ref. [1].) If the /-functions are chosen according to eq. (2.8), the 

y the Hankel det 

In • • • fn+N-1 



T function t„ is then given by the Hankel determinant 






(2.10) 



fn+N-l ■ ■ ■ .fn+2N-2 

for 1 < A^ < M. It should be noted that, even when the set of functions {fn }fLi is chosen 
according to Eq. (2.8), no derivatives appear in the r function, and therefore Eq. (2.10) cannot 
be considered a Wronskian in the same sense as for the KP equation (cf. Eq. (1.9) in Ref. [1]). 
Nonetheless, this choice produces a similar outcome as in the KP equation. Indeed, similarly 
to Ref. [1], we have the following: 

Lemma 2.1. Let /„ be given by Eq. (2.9), with 0n (j ~ I,. . . , M) given by Eq. (2.7). Then, 
for 1 < N < M — 1 the t function defined by the Hankel determinant (2.10) has the form 

tW= Y, A(^l,...,^^,) expl'f^ei^^)], (2.11) 

l<il<---<ijv<M \j = l 

where A(ii, . . . , i^v) is the square of the van der Monde determinant, 

A(n,...,ZAr) == ]J (pi^-pi^f. 
l<j<l<N 

Proof. Apply the Binet-Cauchy theorem to Eq. (2. 10), as in Ref. [1]. D 

An immediate consequence of Lemma 2. 1 is that the r function tA is positive definite, 
and therefore all the solutions generated by it are non-singular. Like its analogue in the KP 
equation [1], the above t function produces soliton solutions of resonant type with web 
structure. More precisely, in the next section we show that, like its analogue in the KP 
equation, the above t function produces an (A^_)_, iV_) -soliton solution, that is, a solution 
with N^ = M — N asymptotic line solitons as n ^ oo and N_ = N asymptotic line solitons 
as n ^ — cx). 

Before we turn our attention to resonant solutions, however, it is useful to take a look at 
the one-soliton solution of the 2DTL equation. Let us introduce the function 

d 

Wn{x,t) ^ —log Tn{x,t) , (2.12) 

SO that the solution of the 2DTL equation is given by 

Vn{x,t) ^ —Wn{x,t) . (2.13) 

ot 
If T„ = e^" + e^" , with 9n given by Eq. (2.7) and pi < p2, then w„ is given by 



kiPi +P2) + Hpi -P2)tanhi(0i'^ - #) 



Pi as X —>■ oo, 
P2 as X -^ — cxo, 

which leads to the one-soliton solution of the 2DTL equation: 

Vn = -liPi P2) (f - ^) sech^il^) - 0(2)) . (2.14) 
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In the x-n plane, this solution describes a plane wave u — $(k • x — w t) with x — {x, n), 
having wavenumber vector k — {k^ , fc„) and frequency uj given by: 

k=(pi-p2, logPl -logP2) =: ki,2, ^= =■■ 1^1.2- 

Pi P2 

The soliton parameters (k, oj) satisfy the dispersion relation ujkx + 2 cosh fc„ — 2 = 0. The 
above one-soliton solution (2.14) is referred to as a line soliton, since in the x-n plane it is 
localized around the (contour) line 9n = On ■ Since in this paper we are interested in the 
pattern of soliton solutions in the x-n plane, we will refer to c = dx/dn as the velocity of the 
line soliton in the x direction. (That is, c = indicates the direction of the positive n-axis.) 
For the soliton solution in Eq. (2.14), this velocity is ci^2, where 

Czj = -(logpj - logpj)/(K - pj) . (2.15) 

3. Resonance and web structure in the two-dimensional Toda lattice equation 

We first consider {N^ , l)-soliton solutions, i.e., solutions obtained when A^ = 1. In particular, 
we start with (2,l)-soliton solutions (i.e., A^ = 1 and M — 3), whose r-function is given by 

0(1) 0(2) g(3) 

Tn = e " +e" +e" . 

with pi < P2 < P3 without loss of generality. The corresponding function Wn describes the 
confluence of two shocks: two shocks for n ^ oo (each corresponding to a line soliton for 
Vn) with velocities ci.2 and C2,3 merge into a single shock for n -^ — oo with velocity ci.3, 
with a J given by Eq. (2.15) in all cases. This Y-shape interaction represents a resonance 
of three line solitons. The resonance conditions for three solitons with wavenumber vectors 
{ki.j I 1 < i < J < 3} and frequencies {uJij | 1 < « < J < 3} are given by 

ki,2 +k2,3 = ki,3, and uji^2 + (^2,3 = ioi^s, (3.1) 

which are trivially satisfied by those line solitons. We should point out that this solution is 
also the resonant case of the ordinary 2-soliton solution of the 2DTL equation. As mentioned 

earlier, ordinary 2-soliton solutions are given by the N ^ 2 r-function (2.4) with (2.6). The 

(2) 
explicit form of the tA -function is 



T, 



(2 



(2,3) 



^ = (Pi-P3)e" +(P1-P4)e" +(P2-P3)e" +(p2-P4)e 



where for brevity we introduced the notation ^n = ^n + 6'„ , and where 0„ is given by 

(2) 
Eq. (2.7), as before. Note that \ip2 = P3, the t„ -function can be written as 

T^2) ^ ge„ +e„ +e„ ^^^^ _ ^^^ ^ ^ e„ + (p^ _ p^) e »- + (p^ -Pi)e "^ ] , 

where A — exp((?g — ^0 ) = constant. Since the exponential factor e^" +^" +^" gives 

(2\ (2) 

zero contribution to the solution Vn — dtdx^ogT„ , the above rA -function is equivalent 
to a (2,l)-soliton solution except for the signs of the phases (more precisely, it is a (1,2)- 
soliton). Note also that the condition p2 = Ps is nothing else but the resonance condition, and 
it describes the limiting case of an infinite phase shift in the ordinary 2-soliton solution, where 
the phase shift between the solitons as n ^ ±00 is given by 

S = {pi- P3){P2 - Pa)/[{P2 - P3){pi - Pa)] ■ (3.2) 

The resonance process for the (A^+, 1) -soli ton solutions of the 2DTL equation can be 
expressed as a generalization of the confluence of shocks discussed above (cf. Ref. [1]). 
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We next consider more general {N^, Af_)-soliton solutions. Following Ref. [1], we can 
describe the asymptotic pattern of the solution in the general case iV 7^ 1 by introducing a 
local coordinate frame (^, n) in order to study the asymptotics for large \n\, with 



X = en + S, ■ 

1 Eq. (2.7) tl 
C^=P^e + '?^(c)n + ^^'^ for i = l,...,M, 



The phase functions 0n in / in Eq. (2.7) then become 



with 

1^^{c) := p,(c + (1/p,) logpj) . 

Without loss of generality, we assume an ordering for the parameters {pi \ i = 1, . . . , M}: 
< pi < p2 < ■ • • < Pm ■ Then one can easily show that the lines r] ~ r/i (c) are in general 
position; that is, each line rj ~ rji (c) intersects with all other lines at M — 1 distinct points in 
the c-T] plane; in other words, only two lines meets at each intersection point. 

The goal is now to find the dominant exponential terms in the t„ -function (2.1 1) for 
n -^ ±00 as a function of the velocity c. First note that if only one exponential is dominant, 
then Wn = dx^ogTn is just a constant, and therefore the solution Vn = dtWn is zero. 
Then, nontrivial contributions to Vn arise when one can find two exponential terms which 
dominate over the others. Note that because the intersections of the rji's are always pairwise, 
three or more terms cannot make a dominant balance for large |ri|. In the case of (A^+, 1)- 
soliton solutions, it is easy to see that at each c the dominant exponential term for n ^ — cxo 
is provided by only 771 and/or tjm, and therefore there is only one shock (N^ = 1) moving 
with velocity cim corresponding to the intersection point of 771 and r/M- On the other hand, 
as n ^ 00, each term rjj can become dominant for some c, and at each intersection point 
rjj = rij+i the two exponential terms corresponding to rjj and rjj+i give a dominant balance; 
therefore there are N+ — M — 1 shocks moving with velocities Cj_j+i for j = 1, . . . , M — 1. 

In the general case, N ^ l,theTri -function in (2. 11) involves exponential terms having 
combinations of phases. In this case the exponential terms that make a dominant balance can 
be found using the same methods as in Ref. [1]. Let us first define the level of intersection of 
the ?7i(c). Note that q^ in Eq. (2.15) identifies the intersection point of 77^(0) and rjj{c), i.e., 
Vt{ct,j) =%(Q,i)- 

Definition 3.1. We define the level of intersection, denoted by aij, as the number of other 
rji 's that are smaller than rji {ci,j ) — rjj (cij ) at c — aj. That is, 

0-»J — \{Vl I VliCt^j) < ?7^(C^J) = ??j(Qj)}| • 

We also define I{s) as the set of pairs (rji, rjj) having the level cjij — s, namely 
I{s) := {(?7j, -qj) I a^^j = s, for Kj]. 

The level of intersection lies in the range < cr^.j < M — 2. Note also that the total number 
of pairs [rji, rjj) is 

- s _. M-2 

f) =iM(M -1)= 5: |/(.)|. 

^ n=0 

One can show that: 

Lemma 3.2. The set I{s) is given by 

I{s) = {(?7i, 77M-s+i-i) I i = 1, . . . , s + 1} . 
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Proof. From the assumption qi < q2 < ■ ■ ■ < qm, we have the following inequality at 
c = Cij (i.e. T]i = rij) for i < j, 

?7i, . . . ,?7i-i, ?7j+i, ...,r]M<Vi'^Vj< Vi+i, ■ ■ -^Vj-i ■ 
Then taking j ~ M — s — 1 leads to the assertion of the Lemma. D 

Now define N- ~ N and iV+ = M — N. The above lemma indicates that, for each 
intersecting pair (rji, rjj) with the level A^^^ — 1 (iV_ — 1), there are N- — 1 terms 77/ 's which 
are larger (smaller) than rji = rjj . Then the sum of those iV_ — 1 terms with either rji or rjj 
provides two dominant exponents in the r„ -function for n ^ 00 (n ^ —00) (see more 
detail in the proof of Theorem 3.3). Note also that |/(A^± — 1)| = A^^. Now we can state our 
main theorem: 

Theorem 3.3. Let Wn be defined by Eq. (2.12), withTn given by Eq. (2.11). Then Wn has 
the following asymptotics for n -^ ±00 : 

(i) For n ^ 00 and x ~ Ci^N_^i n + ^ for i ~ 1, • . • , N^ , 



^.(+,+):=Ef=t;iP. asC 



Wn ' % ,, . ^ ^^"AT+i-l 



(ii) For n —>■ —00 and x = Ci^N++i ^ + C foi" i = 1, ■ ■ ■ , -/V_ , 



— 00 . 



\ K,{+, -) := X;]=l Pj + J2f=i^^ PM-j+1 as C ^ 00 , 

\ Ki{-, -) := J2]=i P'j + Y.j^i PM^3+t as C ^ -00 . 
where Cij is given by Eq. (2.15). 

Proof. First note that at the point rji = rjM+i, i.e., (77^, rj^^i) E /(iV+ — 1), from Lemma 3.2 
we have the inequality, 

Vi = VN+i < '7i+1i'7j+2, • • • ,rii+N-i ■ 



N-1 

This implies that, for c — — (logpTv+i — ^oSPi)/{PN+i — Pi), the following two exponential 
terms in the t„ -function in Lemma 2.1, 

N+i-l s / N+i 



exp (^ Y. eiPj , exp (^ ^ ^(f) j 



j=i ^ ^ j=i+l 

provide the dominant terms for n ^ 00. Note that the condition rji ~ rjM+i leads to 
c = Ci^N+i = —{iogPN+i — ^o&Pi)/{PN+i —Pi)- Thus the function w„ can be approximated 
by the following form along x ~ Ci^^+i n + £_ for n ^ 00; 

Wn - ^log(A,(-,+)e^^(-+)? + A,(+,+)e^^(+'+)«) 

A,(-, +)e^.(-^-)« + A,(+, +)ei<d+.+K 
^ K,{-, +)A,(-, +)e(P'-P"+')^ + K,i+, +)A,(+, +) 
A,(-, +)e(P'-P"+-)f + A,(+, +) 
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where 

A,(-,+) = A(j,...,7V + z-l) cxpf J2 ^o'M 

A,(+,+) = A(* + l,...,iV + j) expf Y. ^o' 

Now, from pi < pN+i it is obvious that w„ has the desired asymptotics as ^ ^ ±00 for 
n ^- 00. 

Similarly, for the case of {rji, rj^^+i) E I{N^ — 1) we have the inequality 

vi7V2,- ■ ■ ,111-1, 'nN++i+i ■ ■ ■ ,riM <'ni = VN++i ■ 

N~l 

Then the dominant terms in the tA -function on x = c^ Ar^+in + ^ for n -^ —00 are given 
by the exponential terms 

E ^i^' + E e^"^'^^^ ) , exp ( E ^i^^ + E e-+^^ ) • 
j=i j=i / \j=i j=i / 

Then, following the same argument as before, we obtain the desired asymptotics as ^ ^ ±00 
for n — > — cxD. 

For other values of c, that is for c j^ Ci^^^ +i and c 7^ Cj jv^+i, just one exponential term 
is dominant, and thus «;„ approaches a constant as \n\ -^ 00. This completes the proof. D 

Theorem 3.3 determines the complete structure of asymptotic patterns of the solutions 
Vn{x, t) given by (2.10). Indeed, Theorem 3.3 can be summarized as follows: As n ^ 00, 
the function w„ has N+ jumps, moving with velocities cj^m^+j for j — 1, . . . , Af+; as 
n -^ —00, Wn has N^ jumps, moving with velocities Q^jv^+i for i ~ 1, • • • , N^. Since 
each jump represents a line soliton for Vnix^t), the whole solution therefore represents an 
(A^+, A^-) -soli ton. The velocity of each of the asymptotic line solitons in the (A^+,iV_)- 
soliton is determined from the c-rj graph of the levels of intersections. As an example, in 
Fig. 1 we show a (2,l)-soliton solution (also called a Y-shape solution, or a Y-junction), a 
(2,2)-soliton solution, a (2,3)-soliton solution and a (3,3)-soliton solution. 

Note that, given a set of M phases (as determined by the parameters pi for i = 
1, . . . , M), the same graph can be used for any (A^+, A^_)-soliton with N^ + A^+ = M. 
In particular, if M := 2A^, we have N+ = N^ = N, and Theorem 3.3 implies that the 
velocities of the A^ incoming solitons are equal to those of the A^ outgoing solitons. In 
the case of the ordinary multi-soliton solution of the 2DTL equation, the r-function (2.4) 
does not contain all the possible combinations of phases, and therefore the theorem should 
be modified. However, the key idea of using the levels of intersection for the asymptotic 
analysis is still applicable. In fact, by considering the t„ -function given by the Casorati 
determinant (2.4) with fi = e^" + e^" for i = 1, . . . ,N and pi < p2 < ■ ■ ■ < P2N, 
one can find the asymptotic velocities for the ordinary A^-soliton solutions as C2i-i,2i = 
— {\0gp2i — \ogp2i-i)/{p2i ~P2i-i) Note that these velocities are different from those of the 
resonant A^-soliton solutions. 

Note also that even when A^+ = N^ — N, the interaction pattern of resonant soliton 
solutions differs from that of ordinary A^-soliton solutions. As seen from Fig. 1, the resonant 
solutions of the 2DTL obtained from Eq. (2.10) are very similar to the solitons of the KP 
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Figure 1. Resonant solutions of the two-dimensional Toda lattice: (a) (2,l)-soliton solution 
(i.e., a Y-junction) at t = 0, with AT = 1, M = 3, pi = 1/4, p2 = 1/2, P3 = 2; (b) (2,2)- 
soliton solution at t = 14, with N = 2, M = A, pi = 1/8, P2 = 1/2, ps = 1 and p4 = 4; 
(c) (3,2)-soliton solution at t = 10, with N = 2, M = 5, pi = 1/10, p2 = 1/5, P3 = 1/2, 
P4 = 1 and p5 = 6; (d) (3,3)-soliton solution at t = 14, with N = 3,M = 6,pi = 1/10, 
P2 = 1/4, p3 = 1/2 and p4 = 1, ps = 2, pg = 4. In all cases the horizontal axis is n and 
the vertical axis is x, and each figure is a plot of q{n, x,t) in logarithmic grayscale. Note that 
the values of n in the horizontal axis are discrete. 



and coupled KP equation [1, 4, 5], where such solutions were called "spider-web" solitons. 
(In contrast, an ordinary iV-soliton solution produces a simple pattern of N intersecting 
lines.) The web structure manifests itself in the number of bounded regions, the number of 
vertices and the number of intermediate solitons, which are respectively (N_ — l)(iV+ — 1), 
27V_7V+ -- M and 3iV_iV+ - 2M for an (iV+, Ar_) -soliton solution [1]. (In contrast, 
an ordinary A^-soliton solution has (N — 1)(N — 2)/2 bounded regions and N{N — l)/2 
interaction vertices.) Finally, it should be noted that, as in the KP equation, only the Y-shape 
solution is a traveling wave solution. All other resonant solutions (as well as ordinary iV- 
soliton solutions with A^ > 3) have a time-dependent shape, as shown in Ref. [1]. 
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4. The fully discrete 2D Toda lattice equation 

The 2DTL equation (2. 1) is a differential-difference evolution equation, since only one of the 
independent variables is discrete, while the other two are continuous. Hereafter, we refer to 
Eq. (2.1) as a semi-continuous case. We now consider a fully discrete analogue of the 2DTL 
equation (2.1), namely 

^l ^niQl.m.n = Vtm-l,n+l — Vl+l^m-l,n — Vl^m.n + Vj+l,m,n-l , (4.1) 
VLm,n = (SKy^ l0g[l + Sk (cxpQ,^„^„ - 1)] , 

with l,m,n ^ Z, I and m being the discrete analogues of the time t and space x coordinates, 
respectively, and where A'^ and A^ are the forward and backward difference operators 
defined by 

// + l,m,n Jl^va.n 






Jl,'m.7i Jl.7n—l,n 



(4.2) 



(4.3) 



Equation (4. 1), which is the discrete analogue of Eq. (2. 1), can be written in bihnear form [2] 
in a manner similar to Eq. (2.2): 

~~ 'H.rn— l,n+l'n+l,7n,n — 1 '^l-\-l,ni—l,n'^l.m,n-j \^-^) 

with Qi^jn^n related to r/^„,„ by the transformation Vi^m,n = A,+ A„ logT;_„,„, i.e., 

(4.5) 



Cl.m 



1 Ti+irn+l,n-lTLm.n+l 
n = log '— 



Note that Qi^m.n = log[l + (e'^'^^'™" - 1)/5k\. Special solutions of Eq. (4.4) (which is the 
discrete analogue of Eq. (2.2)) are obtained when the r function Ti^m.n is expressed in terms 



AN) 



of a Casorati determinant Ti^m.n — Timn ^^ [^1 

^ Ji. 

f(2) 
AN) 



J l.m^n J l.m,n-{-l 
„(2) „(2) 

J l,m,n -z /,m,n+l 



/•(I) 

Jl m.n+N-l 



fl 



l,m,n+N~l 



AN) AN) 

'' l.m,n •> l,m.n-{-l 



(4.6) 



1,2,- 



where each of the functions {/;„„, 
dispersion relations: 

If we take as a solution for Eqs. (4.7) and (4.8) the functions 



AN) 

Jl,m,n+N-1 

, N} satisfies the following discrete 



// 



(^) 



<P{pt) + (f>{qt 



with 



the T function (4.6) yields a iV-soliton solution for the discrete 2DTL Eq. (4.4). 



(4.7) 
(4.8) 

(4.9) 
(4.10) 
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As in the semi-continuous 2DTL, however, solutions of Eq. (4.4) can also be obtained 
when we consider the r function defined by the Hankel determinant 

Jl,m^n //,m,n+l ' ' ' Jl.7n,n+N —1 

Jl,m.n-\-l //,?n,n+2 * ' * Jl^va.n+N 



SN) 



Jl,m,n+N-1 Jl,m,n+N 



fl,m,n+2N- 



(4.11) 



where 



M 



(4.12) 



(4.13) 



fl,fn,n = ^aii(f){pi) , 

which corresponds to choosing 

for i = 1,. . .,N. Without loss of generality, we can label the parameters pi so that 
< pi < P2 < ■ ■ ■ < Pm-1 < Pm- Then, as in the semi-continuous 2DTL, we have 
the following: 

Lemma 4.1. Let fi.m,n be given by Eq. (4.13). Then, for 1 < N < M — 1, the t function 
defined by the Hankel determinant (4.11) has the form 

N 

rW= J2 A(^l,...,^^^)^a.,<^(K,) 

l<n<---<ijv<M j=l 

where A{ii, . . . ,1^) is the square of the van der Monde determinant, 

A{ii,...,iN) = Jl (pi^ -; 
i<j<;<jv 



-po^. 



(4.14) 



(4.15) 



Proof. Again, the result follows by applying the Binet-Cauchy theorem to the Hankel 
determinant (4.11). D 

Unlike its counterpart in the semi-continuous 2DTL equation, the t function in Eq. (4. 1 1) 
cannot be written in terms of a Wronskian, since no derivatives appear. However, as in the 
semi-continuous 2DTL equation, the t function thus defined is positive definite, and therefore 
all the solutions generated by it are non-singular. In the next section we show that, like 
its analogue in the semi-continuous 2DTL equation, the above r function produces soliton 
solutions of resonant type with web structure, and we conjecture that, like in the continous 
case, an (iV+, Af_) -soliton with N+ ^ M - N and N^ ^ N in created. 

Like with the semi-continuous 2DTL equation, however, before discussing resonant 
solutions it is convenient to first look at one-soliton solutions of the fully discrete 2DTL 
equation. Let us introduce the analogue of Eq. (2.12), namely the function 

'n,?n+l,n — 1 



Wl^ra,n = log ■ 



'^l,m,n 



(4.16) 



so that the solution of the discrete 2DTL equation is given by 



l,rn,n = log = W; + l,m,n " Wl,m,n+1 ■ (4.17) 



'^l+l ,711 .n^l .m+l .n 



It is also useful to rewrite the function 



in Eqs. (4.10), (4.12) as 4>{pi) = e*''.™.", where 



l.m.n 



nlogpi - mlog(l + Kp- ^) + /log(l + Spi) + 6*^*^ , 



(4.18) 
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If r„ = e '■"•■" + e i.™.", with pi < p2, then wi^m.n is given by 

— logpi — log(l + npi ) as n — > cx), 

— logp2 ^ log(l + np2 ) as n ^ — oo, 

which leads to the one-soHton solution of the discrete 2DTL equation. In the n-m plane, this 
solution describes a plane wave cx^{Qi.m,n) — *&(k • x — cj/) with x = {n,m), having 
wavenumber vector k = (fc„, km) and frequency uj given by: 

k = (logpi - logp2, - log(l + Kpr^) + log(l + KP^^)) =: ki,2, 

U = - log(l + Spi) + log(l + 6p2) =■■ Wi,2- 

The soliton parameters (k, oj) now satisfy the discrete dispersion relation (e^'^ — 1)(1 — 
e^'^'") ~ (5^(6^*='"+'^" ^g-w-fe„ _ g-w-fe™ _ 2")^ The one-soliton solution (4.17) is referred 
to as a line soliton since, like its semi-continuous analogue, it is localized around the (contour) 
line 0\l^^ — 9\ ^ „ in the n-m plane. Again, we will refer to c — dn/dm as the velocity of 
the line soliton in the n direction. For the above line soliton solution, this velocity is given by 
Ci,2, where now 

C^j = (log(l + np-^) - log(l + Kp-i))/(logp, - \ogpj) . (4.19) 

5. Resonance and web structure in the discrete 2D Toda lattice equation 

As in the semi-continuous case, we first consider (A^-, l)-soliton solutions, i.e., solutions 
obtained in the case A^ = 1, and in particular we start from (2,l)-soliton solutions (i.e., the 
case A^ = 1 and M = 3), whose r-function is given by 

/,(!) /,(2) fl(3) 

with 6*; ^ ^ given by Eq. (4.18), and where pi < P2 < P3 without loss of generality. As 
in the continous case, this solution describes the confluence of two shocks: two shocks for 
m —* oo (each corresponding to a line soliton for Qi^m.n) with velocities ci.2 and ci.3 merge 
into a single shock for m -^ — cxo with velocity ci.3, where c^ j is given by Eq. (4.19) in all 
cases. This Y-shape interaction represents a resonance of three line solitons. The resonance 
conditions for three solitons with the wavenumber vectors {k^ j | 1 < « < J < 3} and the 
frequencies {oji^j | 1 < « < i < 3} are still given by Eq. (3.1), and again are trivially satisfied 
by those line solitons. Furthermore, this solution is also the resonant case of the ordinary 
2-soliton solution of the discrete 2DTL equation, arising in the limit of an infinite phase shift. 
The resonance process for the {N^, 1) -soliton solutions of the discrete 2DTL equation can be 
expressed as a generalization of the confluence of shocks discussed earlier. 

Next we consider more general (A^+, A^_)-soliton solutions. Following Ref. [1] and the 
semi-continuous case, we can describe the asymptotic pattern of the solution by introducing 
a local coordinate frame (^, m) in order to study the asymptotics for large \m\ with 

n = cm + ^ . 
Then the phase functions Ofl^ ^ become 

&'i%,n^i^ogp, + r^,{c)m + 9^'\ for 2 = 1,...,M, 
with 

?7j(c) := logpi (c - log(l + Kp^'^)/\ogpt) . 
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Without loss of generality, we assume an ordering for the parameters {pi \ i — 1, . . . , M}: 
< Pi < P2 < ■ ■ ■ < Pm ■ Then, as in the semi-continuous case, one can easily show that 
the lines rj — rji (c) are in general position. As before, the goal is then to find the dominant 
exponential terms in the T; ^'^^^-function (4.14) for m -^ ±oo as a function of the velocity c. 

First note that if only one exponential is dominant, then wi_m,n — ^^&{'''im+i n-i/umn) 
is just a constant, and therefore the solution Qi^m.n — W(+i^m.n — wi.m.n+i is zero. Then, 
as in the semi-continuous case, nontrivial contributions to Qi^m.n arise when one can find 
two exponential terms which dominate over the others. Also, as in the semi-continuous case, 
since the intersections of the rji's are always pairwise, three or more terms cannot make a 
dominant balance for large \m\. For {N^, l)-soliton solutions, it is easy to see that at each c 
the dominant exponential term for m -^ — oo is provided by only r/i and/or tjm, and therefore 
there is only one shock (N = 1) moving with velocity cim corresponding to the intersection 
point of rji and r/M- On the other hand, as m ^ oo, each term rjj can become dominant for 
some c, and at each intersection point rjj = ryj+i the two exponential terms corresponding to 
rij and rjj+i give a dominant balance; therefore there are A^+ ~ M — 1 shocks moving with 
velocities c^.j+i for j = 1, . . . , M — 1. 

In the general case, N =^ 1, the r^^^^^ -function in (4.14) involves exponential terms 
having combinations of phases, and two exponential terms that make a dominant balance 
can be found in a similar way as in the semi-continuous case. We define again the level of 
intersection of the rji (c). Again, Cij in Eq. (4. 19) identifies the intersection point of r/i (c) and 
T]j{c), i.e., ?7i(cij) = Vjicij)- 

Definition 5.1. We define the level of intersection, denoted by aij, as the number of other 
rji 's that at c ^ Ci_j are smaller than rji (cij ) = rjj {ci_j ). That is, 

o'ij := \{'ni I mici.j) < Vi{ci,]) = Vj{ci,j)}\ ■ 
We also define I{s) as the set of pairs (rji, rjj) having the level (7i_j ~ s, namely 
I{s) := {{■r]^, 7]j) I o-jj = s, for i <j}. 

As in the semi-continuous case, one can then show the following: 
Lemma 5.2. The set I{s) is given by 

^(s) = {{Vi,VM-s+t-i) I J = l,...,s + l}. 

Proof. See the proof of Lemma 3.2. D 

As in the semi-continuous case, let iV_ ~ N and N+ ~ M — N. The above lemma indicates 
that, for each intersecting pair [rji^rij) with the level N+ — 1 {N — 1), there are iV_ — 1 terms 
rji 's which are larger (smaller) than rji = rjj . Then the sum of those iV_ — 1 terms with either 

rji or rjj provides two dominant exponents in the Tj^^^ -function for ra ^ oo [ra ^ — oo). 
We then have the following: 

Theorem 5.3. Let wi^rn,n be afunctiondefinedby Eq. (4.16), withrj^ ^ ^givenby Eq. (4.14). 
Then wi_m,n has the following asymptotics for m — > ±oo : 

(i) For TO — > oo and n — Ci_M-+i to, + ^ for i — 1, . . . , iV^ , 

, , "'V ■ ' ■-) := Ef=t+i logPj as e ^ oo , 

= L,j=t ^Ogpj as ^ ^ -OO . 
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Figure 2. Resonant solutions of the fully discrete two-dimensional Toda lattice: (a) (2,1)- 
soliton solution (i.e., a Y-junction) atl = 0, with pi = 1/10, p2 = 1/2, P3 = 10, (b) (2,2)- 
soliton solution at I = 40, with pi = 1/10, P2 = 1/2, P3 = 2, p4 = 15; (c) (3,2)-soliton 
solution at I = 80, with pi = 1/20, p2 = 1/2, pa = 2, p4 = 10, ps = 60; (d) (3,3)-soliton 
solution at I = 80, with pi = 1/20, p2 = 1/2, p3 = 2, p4 = 10, ps = 20, pa = 120. In all 
cases 5 = K = 1/4; the horizontal axis is n and the vertical axis is n, and each figure is a plot 
of Ql m n in logarithmic grayscale. Note that the values of both m and n in the horizontal 
and vertical axes are discrete. 



(ii) For m 



wi. 



id n = Ci,N^^i m + ^ fori = 1, • • • , N^ , 
^7V-l - 






j+i 



as ^ ^ CX3 , 
as ^ ^ — oo . 



where Cij is given by Eq. (4.19). 
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Proof. Once the obvious modifications are made, the proof proceeds exactly like in the 
semi-continuous case, namely Theorem 3.3. D 

Like its counterpart in the semi-continuous case. Theorem 5.3 determines the complete 
structure of asymptotic patterns of the solutions Qi^m.n given by (4.1 1). Indeed, Theorem 5.3 
can be summarized as follows: As m ^ oo, the function wi.m,n has N+ jumps, moving with 
velocities cj,n_+j for j ~ 1, • • • , N+; as m ^ — cxo, wi^rn,n has N- jumps, moving with 
velocities Ci,M++i for i = 1, • ■ • , N^. Since each jump represents a line soliton of Qi^m.n, 
the whole solution therefore represents an (iV+, 7V_)-soliton. The velocity of each of the 
asymptotic line solitons in the {N+ , N^ ) -soliton is determined from the c-r] graph of the levels 
of intersections. Note that, given a set of M phases (as determined by the parameters pi for 
i — 1, . . . , M), the same graph can be used for any (A^+, Af_)-soliton with N^ + N+ = M. 
As an example, in Fig. 2 we show a (2,l)-soliton solution, a (2,2)-soliton solution, a (2,3)- 
soliton solution and a (3,3)-soliton solution. 

In particular, if M = 2N, we have N+ ~ N^ — N, and Theorem 5.3 implies that the 
velocities of the N incoming solitons are equal to those of the N outgoing solitons. In the 
case of the ordinary multi-soliton solutions of the discrete 2DTL equation, the r-function (4.6) 
does not contain all the possible combinations of phases, and therefore Theorem 5.3 should be 
modified. However, as in the semi-continuous case, the idea of using the levels of intersection 
is still applicable, and one can find that the asymptotic velocities for the ordinary A^-solitons 
generated by the Casorati determinant (4.6) with fi ~ e''" + e^" i = 1, . . . ,N and 
Pi <P2<--- <P2JVareC2i-l,2i = {log{l + Kp^i^)-\og{l + Kp~^^_^))/{\ogp2i-\ogp2i-i), 
Note that, like in the semi-continuous case, these velocities are different from those of the 
resonant A^-soliton solutions. 

The resonant solutions of the fully discrete 2DTL provide the basis for the construction 
of the resonant solution of the ultra-discrete 2DTL, as is shown in the next two sections. 

6. The ultra-discrete two-dimensional Toda lattice 

We now turn our attention to an ultra-discrete analogue of the 2DTL equation. Using 
Eqs. (4.2) and (4.3), we first write the 2DTL Eq. (4.4) in bilinear form as 

(1 — Sk) TiJ^l^m,nn,m+l,n — Tl+l,m+l,nTl^m,n + Sk Tl^m,n+lTl+l,m+l,n-l = , (6.1) 

We define the difference operator A' as 

A'=e-^"(A+-A+)(A+-A+), (6.2) 

where from here on the symbols A^, A+ and A+ will be used to denote the difference 
operators 

A+-e^'-l, A+=c^'"-l, A+ = c^"-l, (6.3) 

and where the shift operators e'^', e'^™ and e^" are defined by e^"/;,„i,„ — fi^m,n+i etc. 
That is, 

A fl^rriAi = .fl+l,m+l,n-l + .fLm,n+l ^ .fl+l,m,n ^ .fLrn+l,n ■ (6.4) 

Using Eqs. (6.3) and (6.2), we can rewrite Eq. (6.1) as 

(1 - (5k) + (5k CXp [A' log Ti^,n,n] = CXp[A+ A+ log Ti^,n,n] , (6.5) 

which becomes, taking a logarithm and applying A' (assuming Sk ^ 1), 

Sk, 



A' log 



1-Sk 



exp(A'logT;,™,„) 



A+A+A'logT,,™,„. (6.6) 
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We now take an ultra-discrete limit of Eq. (6.6) following Refs. (2.1) [9, 10]. This is 
accomplished by choosing the lattice intervals as 

6e - 0-"-/% K, = e-^/% (6.7) 

where r,s (^ Z>o are some predetermined integer constants, and by defining 

<m,« = A'elogT;';„„, (6.8) 

Taking the limit e ^ 0+ in Eq. (6.6) and noting that lim£_,o+ e log(l + e"^/^) = max(0, X), 
we then obtain 

A^A+ W(,„^„ = A' max(0, «;,„,„ -r - s), (6.9) 

where vi mn = lini vf That is, using Eq. (6.8), 

Wi,m,„ = A' lim £Tf^^. (6.10) 

Equation (6.9) is the ultra-discrete analogue of the 2DTL equation, and can be considered a 
cellular automaton in the sense that vi^m,n takes on integer values. 

Let us briefly discuss ordinary soliton solutions of the ultra-discrete 2DTL Eq. (6.9). As 
shown in Refs. [9, 10], soliton solutions for the ultra-discrete 2DTL Eq. (6.9) are obtained by 
an ultra-discretization of the soliton solution of the discrete 2DTL Eq. (4.4). For example, a 
one-soliton solution for Eq. (4.4) is given by 

T/,m,„ = 1 +??1, (6.11) 

with 

and where 

4>{p) ^ p'\l + 5p)\l + np-^)-"' (6.13) 

as before. We introduce a new dependent variable 

Pf,m,n ==£logTi,m,«, (6.14) 

and new parameters Pi,Qi,Ai e Z as 

e^i/^=Pi, e^^'^-^qi, e^i/"=ai. (6.15) 

Taking the limit e ^ 0+, we then obtain 

Pi,m,n = max(0,ei), (6.16) 

where 

<d^ = A^ + n (Pj - Qi) + ;{max(0, Pi - r) - max(0, Q^ - r)} 

+?7i{max(0, —Qi — s) — max(0, —Pi — s)} , (6.17) 

with e^''/^ = 5 and e""/"^ = k as before, and where pi^m.n — lim£_,o+ Pim n- According to 
Eqs. (6. 10) and (6.14), the one-soliton solution for Eq. (6.9) is then given by 

fLm.n — Pl + l,m+l,n-l + Pl,mAi+l — Pl + l,m,n — Pl,m+l,n ■ (6.18) 

Using a similar procedure we can construct a two-soliton solution. Equation (4.4) admits 
a two-soliton solution given by 

T/,m,„ = 1 +?7l +??2 +^12'?l?72, (6.19) 

with 

(P2 -Pi)iqi -92) 



712 



(gi -P2){q2 -piY 



(6.20) 
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Figure 3. A (2, l)-resonant soliton solution (i.e., a Y-junction) for the ultradiscrete 
2DTL Eq. (6.9), with Pi = -5, P2 = 1, P:i = A, r = 3, s = I, Ai - A3 = 1, 
A2 — As = 1. From left to right, the four plots represent the solution respectively atl = 0, 
I = 2,1 = 4 and I = 6. The dots indicate zero values ofvi^„. 



and where rji = ai(j>(pi) / (l){qi) as before (i — 1, 2). In order to take the ultra-discrete limit of 
the above solution, we suppose without loss of generality that the soliton parameters satisfy 
the inequality 

0<pi <p2<q2<qi- (6.21) 

= elogrj,„i^„, as well as integer parameters 



Introducing again the dependent variable pf. 

Pi, Qi and Ai as 



.^./e 



e'?-/^ = ft, 



e-4-/^ = a,. 



(6.22) 



e "- =pi, 
(i — I, 2), and taking the limit of small e, we obtain 

Pi,m,n = max(0, 61, 62, 61 + 62 + P2 - Q2), (6.23) 

where 6^ (i = 1,2) was defined in Eq. (6.17), with pi^m^n — limg^o+ Pf m n ^g^i"' ^"^ 
where vi,m,n is obtained from pi,m,n using Eq. (6.18). Note that Pi < P2 < Q2 < Qi- 

More in general, starting from Eqs. (4.6) and (4.9) (with < pi < P2 < • • • < Pm < 
Qm < Qn-i < ■ ■ ■ < qi) and repeating the same construction, one obtains the A^-soliton 
solution of the ultra-discrete 2DTL Eq. (6.9) as [10] 



Pl,r. 



max 



J2 M*©i+ J2 PiPi'iP^' -Q^') (6-24) 

l<i<M l<i<i'<M 

where max^=o.i indicates maximization over all possible combinations of the integers /i^ = 
0, 1, with i = 1, . . . , M. Again, vi^m.n is obtained from pi^m.n via Eq. (6.18). 

Ordinary soliton solutions corresponding to the above choices were presented in 
Ref. [9, 10]. In the next section we show how this basic construction can be generalized 
to obtain resonant soliton solutions. 



7. Resonance and web structure in the ultra-discrete 2D Toda lattice equation 

Following Ref. [1], we now construct more general solutions of the ultra-discrete 2DTL 
equation (6.9) which display soliton resonance and web structure. 

We first consider the case of a (2, l)-soliton for Eq. (4.4), which is given by 

r/,™,„=a+C2+6, (7.1) 
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Figure 4. An L-shape (2, l)-resonant soliton solution for Eq. (6.9), with Pi = — 5,P2 = — 1> 
P3 = 4, r = 3, s = 1, Ai - ^3 = 1, ^2 - A3 = 1. From left to right, the four plots 
represent the solution respectively aX I = 0, 1 = 2, 1 = 4 and / = 6. 



where 

^i ^ ai(j){pi) (7.2) 

(i^ 1,2,3), with 

Hp)^p''{l + 5pyil + ^p-Y"' (7.3) 

as before, and where again we take < pi < P2 < Ps- As in the previous section, we 
introduce the new dependent variable 

Pf,m,n =£l0g^/,m,". (7-4) 

and new parameters as 

e^*/^=P„ e-^.A^a, (7.5) 

(i — 1, 2, 3), with e^'/^ = 6 and e^*/^ = k as before. Taking the limit e — > 0+, we then 
obtain 

Pl,Tn,n = max(Ji>i,i?2,-R3) (7-6) 

where pi^m.n — limg^o+ Pf m n ^^ before, but where now 

i?j = Aj + nPi + I max(0, Pi - r) - m max(0, -Pj - s) {1.1) 

(i = 1,2, 3). Figure 3 shows that this solution, which again can be called a (2, l)-soliton, is a 
Y-shape solution. We note that, interestingly, an L-shape solution can be obtained instead of 
Y-shape solution for different solution parameters. An example of such a L-shape soliton is 
shown in Fig. 4. No analogue of this solution exists in the 2DTL and its fully discrete version. 
Next, we consider the case of a (2, 2)-soliton for Eq. (4.4) following Ref. [1]. Let us 
consider the following r function 

Jl.m,n //,m,n+l 

//,?n,n+l Jl,'m.n-\-2 



Tl,^ 



(7.8) 



where 



Am^n = a + 6 + 6 + ^4 , (7.9) 

where ^ (i = 1, . . . , 4) is again defined as in Eq. (7.3), and where < pi < P2 < Ps < Pi 
holds. We introduce again the new parameters c^*/^ — pk and e^''^^ ~ at {k ~ 1, ... ,4) 
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Figure 5. Snapshots illustrating the temporal evolution of a (2, 2)-resonant soliton solution 
for Eq. (6.9), with Pi = -7, P2 = -5, P3 = 1, P4 = 3, r = 4, s = 2, Ai = -10, 
A2 = -6, A3 = 0, A4 = 2: (a) I = -12; (b) I = -7; (c) I = A; (d) I = 7. As in Figs. 3 
and 4, the horizontal axis is n and the vertical axis is m. Since the interaction extends over a 
wider range of values of m and n, the solution is now plotted in grayscale, in a similar way as 
in Figs. 1 and 2; the values of t);_m,n however are still discrete, as in Figs. 3 and 4. 



and the new dependent variable p[ 
obtain 



elogr;.m^„. Taking the limit e -^ 0+, we then 



where pi . 



Pl,m,n 

lim. 



1<2<J<4 



(7.10) 



A'y = Rt 



^0+ Pi.m,n^ ^s before, and 

i?,, (7.11) 

and with Rj given by Eq. (7.7) as before. Figure 5 shows the temporal evolution of a (2, 2)- 
soliton solution. Note the appearance of a hole in Fig. 5. 

Like in the 2D Toda lattice (2.1) and its fully discrete version (4.1), we now consider 
more general resonant solutions for the ultra-discrete 2DTL (6.9). We start from the general 
T function defined in Eq. (4.11), and introduce again the parameters e^''/'^ — pk and 



=A,/e 



ttfc (fc — 1,2,..., M) and the variable p[, 



£logTi_„_„, together with 



e '■/'^ = (5 and c '^/'^ = k. Taking the limit e -^ 0+, we then obtain the following solution of 
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Figure 6. Snapshots illustrating the temporal evolution of a (3, 3)-resonant soliton solution 
for Eq. (6.9), with Pi = -10, P2 = -7, P3 = -5, P4 = -1, P5 = 4, Pe = 5, r = 7, 

s = 4, Ai = -8, A2 = -6, A3 = 0, A4 = 2, As = 4, As = 7: (a) I = -15, (b) 
I = -10, (c) I = 0, (c) I = 10. 



the ultra-discrete 2DTL (6.9): 



max 

l<il<-<iN<Af 



N N 

^i?.^. +2^(j-l)P,^. , (7.12) 

Lj = l J=2 

where again lim£^o+ Pi mn = Pi,m.n, with the maximum being taken among all possible 
combinations of the indices ij {j ^ I, . . . , N), and where once more we have 

i?j = Ai +nP.i + ?max(0,Pi - r) - mmax(0,-Pj - s) . (7.13) 

Equation (7.12) produces complicated soliton solutions displaying resonance and web 
structure. As an example, in Fig. 6 and Fig. 7 we show some snapshots of the time evolution of 
a (3, 3)-resonant soliton solution and a (4, 4)-resonant soliton solution. Indeed, we conjecture 
that, similarly to its counterparts for the 2DTL and in its fully discrete analogue, Eq. (7.12) 
yields the (iV+, iV_)-soliton solution of the ultra-discrete 2DTL Eq. (6.9), with N^ = N 
and Nj^ = M — N . Unlike the semi-continuous and fully discrete cases, however, we were 
unable to prove this conjecture using the techniques introduced in Ref. [1]. In this respect, 
it should be noted that solutions of the ultra-discrete 2DTL arise as a result of the properties 
of the maximum function, and therefore their study might require the use of techniques from 
tropical algebraic geometry, which is a subject of current research [13, 14, 15]. 
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Figure 7. Snapshots illustrating the temporal evolution of a (4, 4)-resonant soliton solution 
for Eq. (6.9), with Pi = -15, P2 = -12, P3 = -9, P4 = -3, P5 = 1, Pe = 1, PV = 4, 

Pg = 7, r = 7, s = 4, Ai = -8, A2 = -6, A3 = 0, A4 = 2, A5 = A, Ag = 7, A7 = 8, 
As = 10: (a) I = -10, (b) i = 0, (c) « = 10, (d) I = 20. 



It should also be noted that the interaction patterns in the ultra-discrete system differ 
somewhat from their analogues in the semi-continuous and fully discrete cases. In particular, 
low-amplitude interaction arms may disappear when considering the ultra-discrete limit. 
Furthermore, the specific interaction patterns in the ultra-discrete limit depend on the value of 
the parameters r and s, and different kinds of solutions may appear for different values of r 
and s. In particular, large values of r and s tend to result in the production of several vertical 
solitons, as shown in Figs. 6 and 7. In order to preserve the soliton count in these cases, 
all the outgoing vertical solitons should be counted as one, as should the incoming ones. In 
this sense, a set of outgoing or incoming vertical lines can be considered as a bound state of 
several solitons. A full characterization of these phenomena and their parameter dependence 
is however outside the scope of this work. 

8. Conclusions 



We have demonstrated the existence of soliton resonance and web structure in discrete 
soliton systems by presenting a class of solutions of the two-dimensional Toda lattice (2DTL) 
equation, its fully discrete analogue and their ultra-discrete limit. Soliton resonance and web 
structure had been previously found for nonlinear partial differential equations such as the KP 
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and cKP systems. Note that the 2DTL is a differential-difference equation, its fully discrete 
version is a difference equation and their ultra-discrete limit is a cellular automaton, therefore, 
our findings show that resonance and web structure phenomena are rather general features of 
two-dimensional integrable systems whose solutions are expressed in determinant form. 

A full characterization of the solutions, including the study of asymptotic amplitudes 
and velocities and the resonance condition was provided both in the semi-continuous and in 
the fully discrete case. Their analogue in the ultra-discrete case, together with an analysis 
of the intermediate patterns of interactions is outside the scope of this work, and remains as 
a problem for further research. Of particular interest is the ultra-discrete 2DTL, where new 
types of solutions such as the L-shape soliton shown in Fig. 4 appear. 

Finally, we note that the class of solutions presented in this work is just one of the 
possible choices that yield resonance and web structure. Just like with the KP and cKP 
equations, the class of soliton solutions of each of the systems we have considered (namely, 
the 2DTL and its fully discrete and its ultra-discrete analogues) is much wider, and includes 
also partially resonant solutions. The solutions described in this work represent the extreme 
case in which all of the interactions among the various solitons are resonant, whereas ordinary 
soliton solutions represent the opposite case where none of the interactions among the solitons 
are resonant. Inbetween these two situations, a number of intermediate cases exist in which 
only some of the interactions are resonant. As in the case of the KP equation, the study of 
these partially resonant solutions remains an open problem. 
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